Forest operations,tree species composition and decline in rainfall explain runoff changes in the Nacimiento experimental catchments,south central Chile

Few long-term studies have explored how intensively managed short rotation forest plantations interact with climate variability. We examine how prolonged severe drought and forest operations affect runoff in 11 experimental catchments on private corporate forest land near Nacimiento in south central Chile over the period 2008–2019. The catchments (7.7–414 ha) contain forest plantations of exotic fast-growing species (Pinus radiata, Eucalyptus spp.) at various stages of growth in a Mediterranean climate (mean long-term annual rainfall = 1381 mm). Since 2010, a drought, unprecedented in recent history, has reduced rainfall at Nacimiento by 20%, relative to the long-term mean. Pre-drought runoff ratios were <0.2 under 8-year-old Eucalyptus; >0.4 under 21-year-old Radiata pine and >0.8 where herbicide treatments had controlled vegetation for 2 years in 38% of the catchment area. Early in the study period, clearcutting of Radiata pine (85%–95% of catchment area) increased streamflow by 150 mm as compared with the year before harvest, while clearcutting and partial cuts of Eucalyptus did not increase streamflow. During 2008–2019, the combination of emerging drought and forestry treatments (replanting with Eucalyptus after clearcutting of Radiata pine and Eucalyptus) reduced streamflow by 400–500 mm, and regeneration of previously herbicide-treated vegetation combined with growth of Eucalyptus plantations reduced streamflow by 1125 mm (87% of mean annual precipitation 2010–2019). These results from one of the most comprehensive forest catchment studies in the world on private industrial forest land indicate that multiple decades of forest management have reduced deep soil moisture reservoirs. This effect has been exacerbated by drought and conversion from Radiata pine to Eucalyptus, apparently largely eliminating subsurface supply to streamflow. The findings reveal tradeoffs between wood production and water supply, provide lessons for adapting forest management to the projected future drier climate in Chile, and underscore the need for continued experimental work in managed forest plantations.     

High-frequency monitoring of hydrological and biogeochemical fluxes in forested catchments of southern Chile

The variability of rainfall-dependent streamflow at catchment scale modulates many ecosystem processes in wet temperate forests. Runoff in small mountain catchments is characterized by a quick response to rainfall pulses which affects biogeochemical fluxes to all downstream systems. In wet-temperate climates, water erosion is the most important natural factor driving downstream soil and nutrient losses from upland ecosystems. Most hydrochemical studies have focused on water flux measurements at hourly scales, along with weekly or monthly samples for water chemistry. Here, we assessed how water and element flows from broad-leaved, evergreen forested catchments in southwestern South America, are influenced by different successional stages, quantifying runoff, sediment transport and nutrient fluxes during hourly rainfall events of different intensities. Hydrograph comparisons among different successional stages indicated that forested catchments differed in their responses to high intensity rainfall, with greater runoff in areas covered by secondary forests (SF), compared to old-growth forest cover (OG) and dense scrub vegetation (CH). Further, throughfall water was greatly nutrient enriched for all forest types. Suspended sediment loads varied between successional stages. SF catchments exported 455 kg of sediments per ha, followed by OG with 91 kg/ha and CH with 14 kg/ha, corresponding to 11 rainfall events measured from December 2013 to April 2014. Total nitrogen (TN) and phosphorus (TP) concentrations in stream water also varied with rainfall intensity. In seven rainfall events sampled during the study period, CH catchments exported less nutrients (46 kg/ha TN and 7 kg/ha TP) than SF catchments (718 kg/ha TN and 107 kg/ha TP), while OG catchments exported intermediate sediment loads (201 kg/ha TN and 23 kg/ha TP). Further, we found significant effects of successional stage attributes (vegetation structure and soil physical properties) and catchment morphometry on runoff and sediment concentrations, and greater nutrients retention in OG and CH catchments. We conclude that in these southern hemisphere, broad-leaved evergreen temperate forests, hydrological processes are driven by multiple interacting phenomena, including climate, vegetation, soils, topography, and disturbance history.     

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Classifying GRB 170817A/GW170817 in a Fermi duration–hardness plane

GRB 170817A, associated with the LIGO-Virgo GW170817 neutron-star merger event, lacks the short duration and hard spectrum of a Short gamma-ray burst (GRB) expected from long-standing classification models. Correctly identifying the class to which this burst belongs requires comparison with other GRBs detected by the Fermi GBM. The aim of our analysis is to classify Fermi GRBs and to test whether or not GRB 170817A belongs—as suggested—to the Short GRB class. The Fermi GBM catalog provides a large database with many measured variables that can be used to explore gamma-ray burst classification. We use statistical techniques to look for clustering in a sample of 1298 gamma-ray bursts described by duration and spectral hardness. Classification of the detected bursts shows that GRB 170817A most likely belongs to the Intermediate, rather than the Short GRB class. We discuss this result in light of theoretical neutron-star merger models and existing GRB classification schemes. It appears that GRB classification schemes may not yet be linked to appropriate theoretical models, and that theoretical models may not yet adequately account for known GRB class properties. We conclude that GRB 170817A may not fit into a simple phenomenological classification scheme.     

Digital photographic photometry of NGC 5253 and NGC 2915

The surface digital photometry for NGC 5253 and NGC 2915 is given. The luminosity profiles and colour gradients are discussed and compared. It is concluded that NGC 2915 differs from NGC 5253, both in the extension and intensity of the star burst event because of the kinematical behaviour of the captured gas.CONICET, Buenos Aires, Visiting Astronomer, Cerro Tololo Inter-American Observatory, La Serena, Chile, operated by NSF     

On a global expansion of the disturbing function in the planar elliptic restricted three-body problem

Starting with a simple Taylor-based expansion of the inverse of the distance between two bodies, we are able to obtain a series expansion of the disturbing function of the three-body problem (planar elliptic case) which is valid for all points of the phase space outside the immediate vicinity of the collision points. In particular, the expansion is valid for very high values of the eccentricity of the perturbed body. Furthermore, in the case of an interior mean-motion resonant configuration, the above-mentioned expression is easily averaged with respect to the synodic period, yielding once again a global expansion of (R) valid for very high eccentricities.Comparisons between these results and the numerically computed exact function are presented for various resonances and values of the eccentricity. Maximum errors are determined in each case and their origin is established. Lastly, we discuss the applicability of the present expansion to practical problems.     

Lunar motion: Theory and observations

Results of several fits of the lunar theory ELP 2000-82B and of Moons' theory of libration are presented. The theories are fitted both to JPL numerical integrations and to LLR observations     

Possible meteoroid streams associated with (69230) Hermes and 2002 SY50

The orbits of (69230) Hermes and 2002 SY₅₀ are similar and the Earth approaches both of them twice: at the end of October the local orbital minimum distances are smaller than 0.007 AU, and at the end of April the distances are smaller than 0.04 AU. This gives us opportunities to observe the meteors associated with these asteroids. Using the geocentric parameters of the orbital close encounters (the theoretical radiants) and our D _N distance function (Valsecchi et al. Mon. Not. R. Astron. Soc. 304 (1999) 743), we searched for meteoroids originated by Hermes and 2002 SY₅₀. A search among 1830 good quality photographic meteors gave negative results: we found no meteor dynamically similar to Hermes or 2002 SY₅₀. In a second search, done in a set of 62150 radio meteors, we applied two methods (M1, M2) and in both cases we found two streams; the streams found with the M1 method had 43 and 30 members, those found with the M2 method had 39 and 14 members. However, these results do not look convincing, due to the small number of common members in the corresponding streams. We therefore conclude that amongst the IAU meteors used in our search there are no compact streams associated with Hermes and 2002 SY₅₀.     

A seasonal climate model of the atmospheres of the giant planets at the voyager encounter time: I. Saturn's stratosphere

A radiative seasonal model which incorporates a multilayer radiative transfer treatment at wave-lengths longward of 7 μm is presented and applied to Saturn's stratosphere. Opacities due to H₂-He, CH₄, C₂H₂, and C₂H₆ are included. Season-dependent insolation is shown to produce a strong hemispheric asymmetry decreasing with depth at the Voyager encounter times, and seasonal amplitudes of 30°K at the poles are predicted in the high stratosphere. The ring-modulated dependence of the insolation and the orbital eccentricity are shown to have a significant effect. Calculations agree closely with the Voyager 1 and 2 radio occultation ingress profiles recorded at 76°S and 36.5°S for CH₄/H₂ = 3.5 + 1.4/? 1.0 × 10^?3;the estimated errors include modeling systematic errors and uncertainties in the occultations profiles. The possible role of aerosols in the stratospheric heating is analyzed. The Voyager 2 egress profile recorded at 31°S cannot be reproduced by calculations. Some constraints on the C₂H₂ and C₂H₆ abundances are derived. The upper portion of the occultation profiles (p < 3mbar) can be matched for C₂H₂/H₂ = 1.0 + 1.3/?0.6 × 10^?7, C₂H₆/H₂ = 1.5 + 1.8/?0.9 × 10^?6 at 76°S and C₂H₂/H₂ = 4 + 6/?4 × 10^?8, C₂H₆/H₂ = 6 + 9/?6 × 10^?7 at 36.5°N. At the northern occultation latitude, the discrepancy with the concentrations derived from analysis of IRIS spectra by R. Courtin, D. Gautier, A. Marten, B. Bézard, and R. Hanel (1984, Astrophys. J.287) can be explained by a sharp variation of the mixing ratios of these gases with altitude in the upper stratosphere. Other interpretations are discussed.